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Abstract 

The paper introduces a modified version of the classical Coupon Collector's Problem entailing exchanges and cooperation 
between multiple players. Results of the development show that, within a proper Markov framework, the complexity of the 
- - - Cooperative Multiplayer Coupon Collectors' Problem can be attacked with an eye to the modeling of resource harvesting and 

sharing within the context of Next Generation Network. The cost of cooperation is computed in terms of exchange protocol burden 
^ . and found to be dependent on only ensemble parameters such as the number of players and the number of coupons but not on 

' the detailed collection statistics. The benefits of cooperation are quantified in terms of reduction of the average number of actions 

before collection completion. 

o '. 

^ \ I. Introduction 

' The classical Coupon Collector's Game is a process in which a player randomly receives coupons corresponding to one out 
of M possible labels and continues playing until she collects a coupon for each possible label (see, e.g., [1]). 

Coupon Collector's Problems (CCPs) deal with the statistical characterization of what happens in a Coupon Collector's 
,_i Game and, even in their simplest form, are the object of many classical and recent investigations. The reason for this interest 
^) is twofold. On one side they offer challenging questions in the field of mathematical statistics (see, e.g., |l2l|l3l||4l||5l@)- 
On the other side they enjoy a plethora of applications mostly derived from general Bayesian capture-recapture method |T| 
specialized to the most diverse fields such as: counting of biological species or phenomena (see, e.g. [8J and |(9J), IP tracking 
. and Internet characterization (see, e.g., |10| and ifTTl ). recognition of the author of a textual documents JT] as well as many 
other information processing/managing tasks (see, e.g., |T2l, |T3l, |T4l, and ifTSl ). 

Despite this large availability of theoretical and practical results, little is present in the Literature about CCPs with more than 
one player (see @| and (|6l for some early treatment). In particular, what is never addressed is the case in which all players 
aim at completing their collection while accepting to exchange coupons with other players. 

Yet, Cooperative Multiplayer CCPs (CM-CCPs) are promising models for the behavior and performance of resource managing 
ly-j mechanisms in nowadays and foreseeable information and communication technology systems, such as resource harvesting 
PsJ ■ and sharing in Next Generation Networks (NGN) but also knowledge dissemination in social infrastructure as well as content 

, distribution in peer-to-peer organizations to name a few. 
y—i • Note, in fact, that NGN will have to cope with the problem of using resources in un-known, un-planned and un-supervised 
' environments. Such an environment can be coped with by means of a resource management that hinges on two key steps: 
resources harvesting and resource negotiation/exchanging. Though resource harvesting can be a trivial, even randomized task. 
J> ] Negotiation and exchanging needs intelligence and cooperation. Intelligence is needed to administer the trade-off between 
K> ' resources acquisition and release, and balance local and global profit, present and perspective gains. Cooperation is needed to 
rS release a resource when this can be reasonably expected to favor global performance. 

Within that general framework, consider the example of the distribution of DRM-protected multimedia content to a group 
of P subscribers. The complete content is made of individually protected and distinguishable pieces, that can be unlocked by 
only one (not pre-determined) final user. The provider pushes the pieces (possibly packed in envelopes containing more than 
one piece) into an anycast routing network that delivers them to subscribers. Due to unknown internal network mechanisms 
users receive random pieces at random times, retain and unlock what they miss to complete the content, and offer other pieces 
to other users through a multicast routing. Other subscribers request and possibly obtain offered pieces by unicast routing. 
Note that the network used by the content provider is not necessarily the same as the one used for content exchange between 
players, thus leading to a noteworthy increase in service deployment flexibility. 

The overall system is made of a content provider, one of more intelligent distribution networks, and a group of active users 
whose behavior can be modeled as a CM-CCP. 

Note that CM-CCPs that model effective and efficient distribution mechanisms must be such that each coupon entering the 
system is either assigned to a player who needs it or discarded if all the players already have it. This is to say either that 
at any instant the number of players that still need a certain coupon is equal to the number of players minus the number of 
times the coupon entered the system, or that the events "the coupon this label enters the system a number of times equal to 
the number of players" and "all the players have the coupon with this label" are equivalent. 

As a final remark, it is sensible to expect that cooperation plays an important role in distinguishing CM-CCP with P > 1 
players from P non-interacting instances of CCP First, as a matter of fact, interaction entails a cost that must be accounted for. 
Second, even limited cooperation yields advantages in the achievement of players' goal. Hence, classical results on single-player 
CCPs would only yield upper bounds on time-to-completion performance figures. 



2 



player player 

a new lot 
enters the game ^J coupon(s) 
retained 

enters the game 



player 



coupon(s) 
transferred 



coupon(s) 
ffered 



coupon(s) 
requested 



t 



[P - 2)4h 

player 




(p - i)-th 
player 



independently 
for each offered 
coupon 



Fig. 1. Protocol for coupon's exchanging. 



II. Problem statement 

To formally define the problem assume that P players collect coupons, each of them labeled with of one out of M different 
labels. The local goal of each player is to acquire one coupon for each possible label. The common goal of the players is to 
ensure that everybody achieves its local goal as soon as possible. 

An activity burst is triggered by the availability of a lot of L coupons having different labels. The incoming lots are drawn 
independently. When a lot enters the system it is assigned to only one player, who is selected independently each time, and 
so that each player has the same probability of being assigned the lot. 

Once a player receives a lot, she retains the coupons she misses and offers the remaining (duplicate) ones to the other 
players. Then, a contention phase ensues in which each offered coupon is randomly assigned to one of the players requesting 
it, each of them having the same probability of being assigned an offered coupon. A final transfer phase closes activities, in 
which assigned coupons are actually transferred between offering and assigned players. All and only requested coupons are 
transferred. No new burst is started once the common goal has been achieved. 

In Figure [U we schematize the steps of the protocol adopted for coupon's exchange. When the coupon fits in the collection of 
the player first receiving it, no action is initiated. On the contrary, when a coupon can be exchanged, the three-phases protocol 
with offer, request and transfer is adopted. Transfers occur after all requests have been collected. 

Other protocols may be devised to regulate coupon exchange between players. For example, though no hint of a true 
cooperative multiplayer structure is present, the mechanism in [4J can be extrapolated to define a priority-based exchange. 

Comparison between different protocol options in terms of the many merit figures that may characterize them (overhead, 
scalability, flexibility, fairness, etc.) is out of the scope of this paper, which does not aim at protocol optimization but at 
introducing cooperation in CCP and analyzing its effects. From this point of view the reasons to choose the simple protocol 
above are at least threefold: 

• it is completely symmetric and thus avoids unfair behaviors a priori; 

• involved entities (the one assigning the lot at the beginning of the burst and the players exchanging coupons) need a very 
small amount of information on the game structure and state: this avoids any significant startup phase (whose cost should 
be accounted for) and easily copes with a varying number of active players; 

• such simplicity does not prevent it from capturing the costs and benefits of cooperation. 

If all players remain in activity until the common goal is achieved even if they have achieved their local goal we indicate it 
as a "continue-on-completion" (coc ) game. If, on the contrary, players that have completed their collections exit the game, we 
indicate it as an "exit-on-completion" (eoc ) game. The two alternative mechanisms are trivially equivalent for P = 1 while 
they are, in principle, different for P > 1. Results will be given for coc games. Relationships between coc and eoc games will 
be developed to extend these results to the latter. 

The most intuitive quantification of the effectiveness of cooperation mechanisms is the average number of bursts needed to 
reach the common goal, that we will indicate with ^gsoai Concentrating on a specific player, another interesting quantity is 
the average number of bursts needed to complete her collection that we will indicate with B***'**. Finally, it may be of interest 
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to know how many bursts are needed, on average, for the first player to complete her collection, a performance figure that we 
will indicate as B^'^^*. 

Beyond this performance figures, three further quantities may be taken into account as cost factors, namely: the average 
number of offered coupons (C°*) that accounts for the traffic due to the offering phase, the average number of requested 
coupons (C'^^'i) that accounts for the traffic due to the response of the players to the offering phase, and the average number 
of transferred coupons (C*'^) that accounts for the traffic actually needed to make cooperation advantageous for players that 
receive coupons they would have not been assigned if they were playing alone. 

III. Costs 

As far as the costs of coc and eoc schemes are concerned, we may prove the following 

Theorem 1: The costs C°^, C^'^'^, and C^^^ of a eoc game are smaller than the same costs of an coc game with the 
same P > I, M, and L. 

Proof: At any given time, no matter whether in a coc or in a eoc game, P' < P players have not yet finished their 
collection. When a new coupon enters the system P" < P' players miss it in their collection and are ready to compete for it, 
if it is offered. 

The probability that it is offered is -Keoc — {P' ^ P")/ P' for an eoc game and tIcoc = {P ^ P")/P for a coc game. From 
P' < P we get TTeoc < T^coc and thus that is smaller for an eoc game. 

Note also that each time a coupon is offered, it triggers a number of requests equal to P". In expectation this implies that 
also C^'°'^ is smaller for an eoc game. 

Finally, each time a coupon is offered, it triggers one transfer if P" > and no transfer if P" = 0. In expectation this 
implies that also C*''' is smaller in an eoc game. ■ 

In the following discussion we will derive expressions for the costs in the coc game that Theorem [T] guarantees to be upper 
bounds for the eoc case. 

Theorem 2: Independently of the statistics of lot drawing, the average number of transferred coupons is 

^tra = M(p_ 1) 

Proof: Up to P coupons of the same type can enter the system and cause a transfer Assume to sample the game at each 
of the corresponding P time instants. 

The j-th time a coupon of a chosen type enters the system, j — 1 players already have it in their collection, hence there 
is a probability (j — i)/P that the coupon is assigned to one of these players that will initiate the exchange and eventually 
produce a transfer 

Therefore, the average number of transfers of coupons of that type is ~ 1)/^ = (-^ ~ 

Since the histories of the M different types of coupon are independent, the total average number of transferred coupons is 
f(P-l). ■ 
Theorem 3: Independently of the statistics of lot drawing, the average number of requested coupons is 

6 ^ ^ 

Proof: As in the proof of Theorem |2] let us sample the game at every entrance of a coupon of a chosen type. 
The j-th time that this happens, j — 1 players already have it in their collection and P — j + 1 do not. 
Hence, there is a probability (j — 1)/P that the coupon is offered to other players thus producing P — j + 1 requests. 
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Since the histories of the M different types of coupon are independent, the total average number of transferred coupons is 
Theorem 4: The average number of offered coupons is 



Therefore, the average number of transfers of coupons of that type is Ej=i(-P - 3 + - 1)/^ = (P^ - l)/6- 
£ histories of the M different types of coupon are independent. 



Proof: Each offered coupon is either transferred (if a player exists missing it in her collection) or discarded (if no player 
needs it). 

The average number of transferred coupons is, from Theorem |2] ^(P — 1). 

As far as the number of discarded coupons is concerned, note that, since B^°'^^ is the average number of bursts needed for 
aU the players to complete their collection, the average number of coupons entering the game is B^°^^L. The total average 
number of discarded coupons is the number of coupons entering the system minus the MP coupons fitting in the collections. 

Putting all together 

= M(p _ 1) + B^°''^L - MP = B^°''^L - ^(P + 1) 

Note that the last expression can be also interpreted as the total average number of coupons entering the game minus the 
average number of coupons received directly by the players needing them, that results to be -^(P + 1). ■ 
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A. Remarks 

The above Theorems highlight some properties of the costs associated to the cooperation mechanism. 

First, C*'''^ and C'^'^ are independent of the probability that a coupon with a specific label enters the game. Their value do 
not change even if different labels appear with different probabilities. Hence, these two quantities are strictly linked to the 
cooperation mechanism per se. 

Actually, this could have been anticipated thinking that requests and transfers depend only on the entrances of coupons that 
are needed by at least one of the players. Since all types of coupons will eventually enter the system, their probabilities cannot 

affect C^'i and C*'''^. 

The same can be said of L, since the player first receiving the lot treats each coupon separately. 

Another feature that is common to C*''* and C'^'^ is the linear dependence on M: doubling the number of labels implies 
doubling these cooperation costs. 

As far as the dependency on P is concerned, note that though the cost due to transfers is proportional to the number of 
players, the effort implied by signalling coupon requests increases quadratically with it since all players needing a coupon 
apply to obtain it but only one of them is selected to receive it. 

Note finally, that by now no conclusive result on is given due to its dependency on B^°'^^ that will be evaluated in the 
next Sections. 



IV. Performance 

In this Section we will assume that all lots are drawn uniformly, i.e., so that each of the possible (^) lots has the same 
probability. 

Note that since they refer to instant at which either only one or all the players have finished, S^""^' and g**""^* of an eoc game 
are equal to the same performance figures of a coc game with the same P, M, and L. 

When L = 1, a much stronger equivalence holds between the statistical characterization of the whole game evolution in the 
coc and eoc case. In fact 

Theorem 5: If L = 1 the probabUity that the coupon entering the system at the beginning of a burst is finally 
retained by a player pre-selected among those that need it in their collection is the same for a coc and an eoc game 
with the same P and M. 

Proof: At any given time, no matter whether in a coc or in a eoc game, P' < P players have not yet finished their 
collection. When a new coupon enters the system P" < P' players miss it in their collection and are ready to compete for it, 
if it is offered. 

Assuming to chose one of these P" players, the probability that she receives the new coupon is 

1 P-P" 1 _ 1 
p p pii pii 

for a coc game, and 

1 P'-P" 1 _ 1 
^p7 ^7 pii pii 

for an eoc game. ■ 

Note that the above strict equivalence cannot hold in general for L > 1. As a counterexample think of a system with P = 3 
players, M = 2 different labels and L = 2 coupons per lot. The average number of bursts needed by the player that is the 
second to complete her collection can be easily computed in the eoc and coc case. 

In a eoc game, the first player that receives a lot completes her collection and exits, hence S**''*** = l. Between the two 
remaining players, the one receiving the second lot is the second to complete her collection. Hence, the number of bursts 
needed to complete the second collection is 2. Only one player remains that completes her collection as soon as the third lot 
enters the system, thus B^°'^^ = 3. 

In a coc game, the first player that receives a lot completes her collection, hence S^""** = 1 as expected. 

The second lot entering the system can be assigned either to one of the two players with no coupon or to the player that 
has completed the collection. In this latter case (that holds with probability ^), the receiving player assigns the two coupons 
independently to the other players with uniform probability. The probability that the coupons are finally assigned to two 
different players is \ + \ = \- Overall, the probability that no further player completes her collection at the second burst is 
W = \ > ^ thus implying that the average number of bursts needed by the player that is the second to complete her collection 
is strictly greater than 2. 

Regardless whether a further player has completed her collection in the second burst, once a third lot comes in the coupons 
are surely distributed to complete all remaining collection and the game finishes thus making gsoai = 3 as expected. 
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A. Markov model 

The aim of this section is to count the number of bursts needed for game completion. Hence, we need to track the game 
evolution from its initial state toward the achievement of the final common goal. 

A straightforward option would be to evolve a vector containing the number of times a certain label enters the system, i.e., 
an Af-tuple of integers taking values from to P and thus assuming one of (P + 1)*^ different values. 

Note that, since all labels are statistically indistinguishable, what really matters is not how many times each individual labels 
appeared during the game. Rather, we may group labels that appeared once, twice, thrice, and so on, and simply count the 
cardinality of each of these subsets. 

Since each of the first P times a coupon with a certain label enters the systems it is assigned to a player that needs it, 
we may characterize the system by recording how many of the M labels are such that a certain number of players has a 
corresponding coupon. 

More formally, we define a integer (P + l)-tuple S ~ {Sp, . . . , So) where Sj is the number of labels for which exactly j 
players have a coupon. Clearly 

S,>0 

so that one of the components contains a redundant information. 

Due to the sum constraint, the number fip{M) of distinct P-tuples corresponding to observable states is equal to the number 
of ways in which the integer M can be decomposed as the sum of P + 1 non-negative integers or, equivalently, the number 
of ways in which one can choose P objects among M + P, i.e., 

a/ + P 



^ip{M) 



p 



Counting with n the number of bursts since the beginning of the game, we have 5(0) — (0, . . . , 0, M). After that, S{n) can 
be any of the integer vectors satisfying (|T). In particular, if S'(n) = (Af, 0, . . . , 0) = goal then the players have reached their 
common goal and the game is over. 

Since the lots are drawn independently of the state currently characterizing the system, the evolution does not depend on 
the past states but on the current one. Hence, the transitions probabiUties 

Ts",S' = Pr{5(n + 1) = S"\S{n) = S'} (2) 

given for any possible S", S" ^ goal are sufficient to describe the overall game evolution. 

Assuming that S{fi) — goal then, the sequence of states S{n) for n = 0,...,fi — lisa non-stationary stochastic process 
whose first-order characterization is given by the probabilities ps'{n) = Pr{S{n) = S'} for each possible S' ^ goal. Using 
the probabilities Q we have that 

Ts",S' PS'{n) 

is the joint probability that S{n + 1) = S" ^ goal and S{n) — S' ^ goal and that 

ps"{n + 1) = ^Ts".s' PS'{n) 

S' 

is the relationship between the first-order characterization of the system after n and after n + 1 bursts. 
As a final remark on our choice of the state characterizing the system note the following 

Theorem 6: If a starting state S' is given, along with an incoming lot, then the system state S" at the end of the 
burst is the same in a coc and in an eoc game. 

Proof: Among the L coupons in the lot, L' < L coupons are missed by at least one player Regardless of the game 
mechanism, those L' coupons will be eventually assigned to one of those players incrementing by one the number of player 
that own them. Since the state S" will take into account the number of players owning each type of coupon, it will be the 
same for coc and eoc games. ■ 

Seen from the point of view of the state S, the dynamic of coc and eoc games is equivalent and so will be any quantity 
computed using only the state evolution. 

The other side of the coin is that, since the state evolution is independent of the completion-exiting policy, the state itself 
cannot entail information about the completion of a strict subset of collections. 
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B. Transition probabilities 

Assume to be in a certain state S{n) ^ goal. The L coupons contained in the lot causing the (n + l)-th burst can be 
partitioned into P + 1 subsets. The /-th subset has cardinahty Ri and contains the coupons for which exactly I players have a 
coupon with the same label. In full analogy with what happens to the state components, we must have < iZ; < min{L, Si{n)\ 
and X)Co = ^■ 

There are (^) equally probable lots. The number of lots in which i?o coupons have a label from So{n) possible labels, 
Ri coupons have a label from 6*1(71) possible labels, and so on, is nfe=o (^/j"^)- Therefore, the probability that the partition 
Ro, Ri, . . . ,Rp applies to an incoming lot is 



-1 p 

n 

fe=0 



/ Sk{n] 

V Rk 



(3) 



Given such a partition it is easy to reason as follows. The Rp coupons in the P-th subset are discarded since no player 
needs them. 

Each of the Rp-k coupons in the (P — fc)-th subset (for k = 1, . . . , P — 1), can be assigned to k players. Given the perfect 
cooperation between players all those coupons will be assigned. 

Hence, since the Rp-k coupons are all different, the number of labels for which P—k+1 players hold a coupon {Sp-k+i{n)) 
increases by Rp-k while the number of labels for which P — k players hold a coupon {Sp-k{n)) decreases by the same 
amount. 

Putting all together, the transitions from S{n) caused by a lot of coupons partitioned into subsets of cardinality Ro,Ri, . . . ,Rp 
leads to an S{n + 1) such that 

f Rp-i if k = P 

Skin + 1) = Skin) +< Rk-i - Rk if < fc < P 
[ -Ro if = 

Conversely, if we know S{n) and Sin + 1) we also know that the lot of coupons causing this transition was partitioned 
into P + 1 subsets with cardinalities 



Rk = 



= < 



J2 Sjin + 1)-Sjir 



P-1 



if 0< A;<P 



ifk = P 



L-^j(5,(n+l)-5,(n)) 

= WiSin + l)-Sin)) 

where the affine function W is defined between (P + l)-dimensional vectors. 

From the expression of W and from the constraints on the Rk we have that feasible transitions are those for which 



0< ^ Sjin + l) - Sjin) <mhi{L,Skin)} 



for < < P and 



(4) 



< L - ^jiSjin + 1) - Sjin)) < min{L, Spin)} 

For this reason, to arrive at a synthetic writing for Ts",s' it is convenient to define the function 

cj>iS",S') = 



= < 



0< Yl S'J -S'^<mm{L,S'„} 



1 if 



j=k+l 



0<L-J2j{Sj -S'^)< mm{L, S'p} 
otherwise 



(5) 
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that evaluates to 1 if the transition from S{n) = S' to S{n + 1) = S" is feasible and to zero otherwise. 
With this, we may recall ^ to write 



fc 



1) Matrix representation: The probabilities Ts".s' of all possible transitions can be arranged into a transition matrix (that, 
with a slight abuse of notation, we will also indicate with T indexed by a pair of integers) by means of a P-dimensional 
embedding of integer tuples into single integers. 

This mapping must disregard one of the component of the state tuple that is redundant. We choose not to consider 5*0. 

This decided, we want to map the set of integer P-tuples {Sp, . . . , Si) such that 5; > and X^^Li — i'^'^o the set of 
integers {0, 1, ... , fip{M) — 1}. 

We indicate with Qp{Sp, . . . , Si) such a mapping in which the function Qp is defined by the following recursion 

Qp{xi, ...,xp)^ 

xi if P = 1 

fJ-p (^^Xk ^ Qp-i{xi, . . . ,Xp_i) if P > 1 

that can be unrolled to give Qp{Sp, . . . ,Si) an explicit expression 

Qp{Sp, . . 




Sp 



Let us prove that Qp is a bijection by induction on P. For P = 1 the fact is trivial. Assume then that P > 1 and that 
Qp-i is a bijection. 

Given an integer g e {0, . . . , iip{M^ — 1} we may compute the xi, . . . , such that q = Qp{xi, . . . , xp) in few steps. 

First, we derive the value of s = J2k=i directly from q — Qp{xi, . . . , xp) = fj,p{s — 1) + Qp^i{xi, . . . ,xp-i). In fact, 
if s were given, the minimum for q would be /ip(s — 1) (obtained for xp = s and xi = for I < P) while its maximum 
would be /ip(s — 1) + /ip-i(s) — 1 (obtained for xi = s and xi = for I > 1). Yet, since by the well-known property of the 
binomial coefficients /ip(s — 1) +/ip-i(s) = fJ'P{s), we conclude that q G {np{s — 1), ... , /ip(s) — 1}. Since for s = 0, 1, . . . 
those ranges are disjoint, s can be directly inferred from q. 

Once that s is known, we may compute Qp_i(a;i, . . . ,xp^i) = q — fJ.p{s ~ 1) and, since Qp-i is a bijection we also 
know xi, .... xp^i. From these we finally get xp — s ~ X^fcLi^ •^fe- 

Beyond being invertible, we also have that Qp{0, ... ,0) = and Qp(goal) = fip{M) — 1 so that it is also a bijection 
between the possible states S ^ goal and the integers {0, . . . ,/ip(A/) — 2}. Hence, we may define the {hp{M) — 1) x 
{fip{M) — 1) transition matrix 

TQj,{S"),Qr{S') = Pr{S{n + 1) = S"\S{n) = S'} 

for all 5', S"'^ goal. 

Note that (|4]i guarantees that, for any feasible transition we have X^jLt+i Sj{n + 1) > Sj{^) for < fc < P so 

that Qp{S{n + 1)) > Qp{S{n)). Since the entries of Tqp(^sii-) Qp(^s'} corresponding to non-feasible transitions are null, the 
matrix T is lower triangular. 

From now on, the same embedding used to arrange transition probabilities in the matrix T will be used also to arrange 
first-order probabihties ps{n) in a one-dimensional array of real numbers also indicated by p and defined as 

PQp(5')H = Pr{^H = 5'} 

In an analogous way, we will adopt the general convention of considering any integer index g e {0, . . . ,/ip(Af) — 2} as 
the representative of a state S — Qp^iq) 7^ goal to allow calculations to be expressed in terms of matrix operations. 

With this, for example, the relationship between the first-order probabilities after n burst and the first-order probabilities 
after n + 1 bursts can be rewritten exploiting standard matrix product as p{n + 1) = Tp{n) — T"p(0). 

To exemplify the construction of the transition matrix T assume first that P = 1. In this case the state is an integer scalar 
simply accounting for the number of coupons in the collection of the unique player and the transition matrix can be written 
straightforwardly. 
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Fig. 2. State embedding and transitions matrix for P = 2, AI = 3, L = 1. 



For P = 1, M = 3 and L = 1 the matrix T is 

/ \ 

r= 1 1/3 

V 2/3 V3 / 

in which the first null row (and thus null eigenvalue) corresponds to the fact that the 0-th state cannot be reached from any 
other state. Note that, as declared before, we are not interested into transitions to the final state that do not appear in the 
matrix. 

If we increase L to 2 the transition matrix reflects this by featuring a further null row corresponding to the new unreachable 
state in which only one coupon entered the game. The matrix is then 

/ \ 

r= 

V 1 V3 V3 J 

When P = 2 the state is the integer pair {S2, Si) whose feasible values can be mapped to the set of integers {0, 1, ... , {M + 
1)(M + 2)/2 — 1}. Since transitions to the final state (Af, 0) are not taken into account, the resulting matrix is an [{M + 
1)(M + 2)/2 - 1] X [(M + 1)(A/ + 2)/2 - 1]. 

For P = 2, M ~ 'i and L ~\, Figure |2}left shows how the two dimensional states different from goalare mapped into the 
integers {0, 1, 2, 3, 4, 5, 6, 7, 8} along with the feasible transitions between the resulting integer states. 

Probabilities for those transitions with the exception of the dashed one are reported in Figure |2}right in the form of the 
transition matrix T. 

For P — 2, M — ?, and L ~ 2, Figure [3}left shows how the two-dimensional states different from goalare mapped into the 
integers {0, 1, 2, 3, 4, 5, 6, 7, 8} along with the feasible transitions between the resulting integer states. 

Probabilities for those transitions with the exception of the dashed ones are reported in Figure [3}right in the form of the 
transition matrix T. 

Note that all the Markov chains we are dealing with are finite and feature only "forward" transition up to a unique absorbing 
state that is the final goal. This class of Markov chains is well-understood and the following results are obtained by specializing 
established methods (see, e.g., [|16i). 
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Fig. 3. State embedding and transitions matrix for P = 2, M = 3, L = 2. 



C. Computation of B^°'^^ 

Theorem 7: If the (/zp(M) - 1) x (/ip(M) - 1) matrix T is such that 

-IP/ 



TQr^{S"),Qp{S') = <I>{S", S') (^^j f[ 

^ ' fc=0 



5;. 



and the numbers r, are defined as 



1 - T,. 



if J = 

■ ^ Tj^uTk if j = 1, . 

/c=0 



.,MP(M) 



then 



oal 



A'p(Af)-2 
j=0 



Proof: If n is such that S(n) = goal then 



B^""*' = E[n] = ^ Pr{n > n} 



n=0 



(6) 



The path is straightforward. At the beginning, no player has a coupon, i.e., the system is in the state 5(0) = (0, 0, . . . , 0, M) 
with probability 1. Since (5p(0, . . . ,0) = we have p(0) = (l ... O) = v, where the column vector v remains 
implicitly defined. 

Starting from this initial condition, the system evolves at each burst and, after n bursts, is characterized by the first-order 
probabilities — T^v. 
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Fig. 4. Normalized performance by varying M for different P and L. 



The probability that n > n is the probability that, after n bursts, the system is in any state but goal, i.e., Pr{n > n} = 



l)p(n) = up{n), where the row vector u remains implicitly defined. 

^goal 



uT"v = u{i - ry^ V 

n=0 



(7) 



where / is the identity matrix. Note that the inversion of (/ — T) is possible since T is lower triangular and thus exhibits its 
eigenvalues on the diagonal, and since such eigenvalues (that are the probabilities that none of the L coupons of the incoming 
lot can be assigned to any player) are less than one if the common goals has not been achieved yet. 

Note finally that / — T is also lower triangular. Hence, the above expression can be easily expanded component-wise to 
realize that, by defining the sequence of r, as in (|6]l, we get the thesis. ■ 

In Figure |4] we report ^ against M for different values of P and L. The adoption of a normalized quantity that 
measures performance in terms of the average number of coupons received directly by each player allows to compare different 
configurations. 

Besides the general and intuitive trend increasing with M, the plots reveal two causes of improvement (i.e., decrease of 

^goal-) 

The first is the availability of lots with L > \. This introduces a correlation between coupon appearance that benefits the 
collection since it reduces the number of duplicates. 

The second is the cooperation between players (i.e., the fact that P > 1) that introduces additional sources of coupons (the 
other P — 1 players) that may contribute to the completion of the collection. 

Note how this second effect is, in this plot, much more evident than the first, to the extent that the case L = 1 is always a 
very good upper bound on f^soai computed for other values of L. 

Actually, this is due to the fact that, for most of the plotted configurations, -g^ ^ 1 thus reducing the effect of coupon 
correlation in lots. 

As a final remark, from ^ we get that the computation of the Tj may greatly benefit from the fact that T is, in general, 
quite sparse. 

In fact, the number of feasible transitions from the generic intermediate state (S'p, S'p_i, . . . , 6*0) is equal to the number of 
ways in which the integer L can be decomposed as the sum of P integers < i?; < min{L, Si). Regardless of the state, this 
number can be easily bounded from above disregarding the constraint Ri < Si. 

Following the same path that led us to compute the total number of feasible state vectors, we obtain that not more than 
(^p^) transitions are feasible from each intermediate state. 

Since there are (^'p^) — 1 states different from the final goal, the number of non-vanishing terms considered in the 



computation of the Tj in ^ is not larger than (^p^) 



1 
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D. Computation of fisoai for P = I 

As noted before, in this case ^i{M) = M + 1 since the state of the system is the degree of filHng of the collection of the 
unique player: from coupons to one coupon for each of M — 1 labels before achieving goal. 
In this case 



j,k 



otherwise 
The matrix T can be written as T = EAE where A is a diagonal matrix with eigenvalues 



A — T ■ 



L / \L 



for j = 0, . . . ,M — I and where 
for j,k = 0,...,AI - 1. 

The matrix E is independent of L and has the noteworthy property of being lower triangular and such that E^^ = E or 

EE = I. 

Based on this, we may compute the terms tj in (|6]l that are the entries of the first column of (/ — T) ^ = E{I — A)^^E, 
i.e. 



i 1 



k=0 ' 

' 'M-k\, ,,M-, (l) f M 



fe=0 

3 



M-k 



Plugging this into (|7]) we get 



ggoal 



M 
L 



v^v^oVV o-ii) 



k=Q \L) \L) j=k 

k=o \lJ \l) 



This expression is a special case of equation (5) in |l3l for equally probable subsets of L out of AI labels 

jjzrj thus reproducing the classical f^soai _ 



Moreover, in the special case L ~ 1, tj turns out to be Tj = -jj^ thus reproducing the classical f^soai _ MJ2fLij ^ 



whose M log M asymptotic trend is well-known. 



E. Computation of B for L = 1 

The computation of this quantity hinges on counting how many coupon assignments are subsumed by each state S and 
distinguishing in how many of them at least one player has finished her collection. This can be done straightforwardly in the 
L = 1 case that is the one addressed here. 

Assume now that n is the number of bursts needed by the first player to complete her collection. As before, the average of 
this quantity can be written starting from its complementary cumulative distribution function as 

oo 

B^""^* ^E[n]^y Pr{n > n} 

In this case Pr{n > n} is the probability that after n burst no player has completed her collection and we may expand 

Pr{n>n}= ^ Pi{n > n\S{n) = S'}Pi{S{n) = S'} 

S'T^^goal 
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Note that 5" spans all the states with the exception of goal since the common goal implies that all players have completed 
their collections. 

To compute Pr{n > n|S'(n) = S'} assume that 5" = {Sp, . . . , So). Corresponding to that state, the X^jLi J^j coupons that 
are in the game, may be assigned to the P players in 

ways, that are all equally probable. 

Among those assignments, there are V in which no player has completed her collection. This number V is the difference 
between U and number of assignments in which exactly j players have completed their collection, for j ~ 1, P. 

Note now that, if 5*0 > then no player may have completed her collection. In general, if 5*^ > then not more than j 
players may have completed their collection. 

Assuming that it is possible, the number of assignments in which exactly j players have completed their collections is the 
number of choices of j players out of P (i.e., (^)) times the number of assignments of the remaining X^fcLi ^Sk — Mj coupons 
to the residual P — j players such that none of them have completed their collections. To count these assignments note that 
since we drop j players that have a coupon for each possible label, if Sk was the number of labels for which fc-players had a 
coupon before, now only k — j out of the remaining P ~ j have a coupon for those labels. 

Hence we may recursively write 

VpiSp,...,So) = Up{Sp,...,So)+ 

max{fc|So — — 0} 



E (^)vp^,{Sp,...,S,) 



3 = 



to yield 



Pr{n > n\S{n) = 5"} - ^^^^'^ 



Up{S') 

Assuming to align all the above conditioned probabilities in the {iip{M) — 1) -dimensional row vector v such that Vj 
Pr{f;, > n\S{n) — Qp^{j)}, a path identical to what was followed for ^Bsoai j^^^j^ 

Relying on the previous ^ we may finally write 



E 

j=0 



F. Computation of B*"^'^ for L = 1 

As before, the derivation hinges on counting how many coupon assignments are subsumed by each state S and distinguishing 
in how many of them the chosen player has finished her collection. This can be done straightforwardly in the L ~ 1 case that 
is the one addressed here. 

Assume now that n is the number of bursts needed by a specified player to complete her collection. As before, the average 
of this new quantity can be written starting from its complementary cumulative distribution function as 

oo 

B*^''' = E[n] = Pr{n > n} 

In this case Pr{n > n} is the probability that after n burst the chosen player has not completed her collection and we may 
expand 

Pr{n >n}^ ^ Pr{n > n\S{n) = S'} Pr{S{n) = S'} 

Note that 5" spans all the states with the exception of goal since the common goal implies that all players have completed 
their collections and that 



Pr{n > n\S{n) = S'} = l-\{ 

k=0 ^ ^ 
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Fig. 5. Global performance compared with performance relevant to eoc games by varying M for different P. 



Assuming to align all the above conditioned probabilities in the {iip{M) — 1) -dimensional row vector ^ such that 
Pr{n > n\S{n) = Qp^ij)}, a path identical to what was followed for ig^ds to 



this 



Relying on the previous ^ we may finally write 



Mp(M)-2 
3=0 



G. Details of performance for L = 1 

In Figure |5] we report , ^-p—, and ^-p— against M for different values of P and for L = 1. As before we measure 
performance in terms of the average number of coupons received directly by each player. 

As expected, we always have B**""^* < 6*'^'^ < S^oai p^^. p = gthis ^ ^first ^ ggoai^ .^j^jjg ^[^pi^ symmeti-y implies 
that B*""'^ = (S**''^* + ^soai) /2 for P = 2. For P = 4 we still have B*""'^ « (B**""^* + ^soai) /2. 

Finally, besides the general and intuitive trend increasing with M, the plots reveal that improvement due to cooperation 
applies to all performance figures. 



V. Conclusions 

This work deals with the statistical characterization of multiplayer coupon collector's games that are a generalization of 
classical coupon collector's games with perspective applications in several Information Technology fields. 

What is addressed is the combination of benefits and costs due to the possibility of cooperation between P players by means 
exchanging of coupons that enter the game in lots each of L different units. The local goal of each player is to complete her 
collection of M distinct coupons. The global goal is the completion of all collections. 

Such an exchange process is regulated by a protocol entailing offer, request and transfer phases. Two playing mechanisms 
are analyzed: one in which players who complete their collection exit the game (eoc games), and one in which they remain 
active and contribute to the overall exchanging (eoc games). 

The average cost of offer (C°*), request (C'°'^) and transfer (C*"^**) phases is computed in analytical terms all yielding very 
simple closed form expressions. The quantities C"''^ and C*'^ turn out to be independent of L and of the statistics of lot 
drawing. 

Costs are larger for coc games than for eoc ones. 
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As far as performance is concerned, an analytical form is given for the average number of activity bursts needed by the first 
player to complete her collection (S**""^*), a chosen player to complete the collection (B^^^^) and by all players to complete 
their collections 

In this case, the equivalence between coc and eoc games is proved as far as B^'^^^ and B^°'^^ are concerned, and when L = 1 
for S**"^. Computations of this merit figures quantifies the effectiveness of cooperation under different points of view. 
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